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Abstract. We have investigated the phenomenon of explosive chromospheric evaporation from an accretion disk 
as a mechanism for fast variability in accreting sources such as low mass X-ray binaries and active galactic nuclei. 
This has been done in the context of advection dominated accretion flows, allowing both high and low states to 
be considered. This mechanism can in principle produce sub-millisecond timescales in binaries and sub-minute 
timescales in active galaxies. However, even considering the possibility that large numbers of these microflares may 
be present simultaneously, the power emitted from these microflares probably amounts to only a small fraction of 
the total X-ray luminosity. 
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1. Introduction 

All models for X-ray binaries (XRBs) involve accretion 
onto compact central objects. Similarly, accretion onto su- 
permassive black holes is the dominant paradigm for the 
engine of active galactic nuclei (AGN). Recently, advec¬ 
tion dominated accretion flows (ADAFs) have been advo¬ 
cated as a natural way of explaining the variety of states 
in which individual binary X-ray sources have been ob¬ 
served (e.g. Chakrabarti & Titarchuk 199f; Esin et al. 
|1997 , 199^ ). Rapid fluctuations have been detected in 
many electromagnetic bands and modeled in a wide vari- 


1995; Tanaka & Shibazaki 1996; Wiita 


ety o f ways (e.g. Krishan fc Wii ta |1994; Wag ner & Witzel 

1996D. In low-mass 


1997, 1998; Tanaka & Shibazaki 1 


Trudolyubov et al. 2001). Seyfert galaxies typically ex¬ 
hibit Lx ~ IQ42-44 gj.g g-i^ with noticeable variations 
observed on t imesc ales as short as min utes (e.g. Lawrence 
& Papadakis 1993 ; Turner et al. 1999 ). 

Many accretion models involve hot low-density coro- 
nae lying over high density, cooler accretion disks (e.g. 
Galeev et al. 11979; White & Holt |1982|; Haardt et al. |1994; 


X-ray binaries, we are presumably dealing with accretion 
onto black holes (BHs) with masses typically 5-15 Mq, 
with luminosities of Lx ~ gj-g g-i |;]^g low/hard 

state, and with fast variations detectable down to 0.01 s 
(e.g. Esin et al. 


Poutanen 1998), and much of the resulting physics is not 


dissimilar to that in ADAF-type models. Recent three- 
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dimensional magnetohydrodynamical simulations of ac¬ 
cretion disks indicate that magnetic fields will be strength¬ 
ened within the disk and that plasma and fields will escape 
from the disk into a strongly magnetized corona (Miller 
& Stone 200C| ). This situation is reminiscent of the con¬ 
ditions prevailing in the solar atmosphere, and therefore, 
our understanding of solar flare related processes may well 
be carried over to accreting sources such as XRBs and 
AGN, which exhibit rapid variations in luminosity (e.g. 
Kuijpers 1995; di Matteo 1998; Krishan et al. 2000). In 
particular, the phenomenon of chromospheric evaporation 
caused by the nonthermal electron heating is often invoked 
to explain the rapid rise of the soft X-ray emission dur¬ 
ing flares and the blueshift of soft X-ray lines of Gaxix 
and Fexxv (e.g. Antonucci et al. 1982). This explosive 
evaporation has been studied using what is known as the 
exploding gasbag model (Fisher [1987 Bastian et al. 1998). 
In this paper, we extend this model to conditions relevant 
to accretion flows in both LMXBs and AGNs, and make 
an approximate evaluation of its possible contribution to 
X-ray emission and variability. 


2. Basic Model 

Figure 1 illustrates the situation we have in mind. We 
consider an ADAF allowing for variations in the accre¬ 
tion rate, M, and thus in the transition radius, rj, which 
is thereby capable of reproducing the spectrum of X-ray 
binaries in both their low/hard state and high/soft states 
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(e.g. Esin et al. 1997 ). The low (high) state is presumed to 
correspond to low (high) values of the accretion rate, nor¬ 
malized to the Eddington accretion rate, m = M/MEdd- 
For larger values of m the ADAF region shrinks (r^ de- 
creases)(e.g. Abramowicz et al. 1996; Narayanet al. 1997D ; 
specific ranges of m are capable of producing the quies¬ 
cent, intermediate and very-high states observed in some 
BH XRBs, such as Nova Muscae 1991 (Esin et al. 19971), 
Cygn us X- 1, GRO J0422-k32 and GRO J1719-24 (Esin 
et al. 199^ , but we shall not be concerned with those de¬ 
tails here. The two temperature plasma within the ADAF 
region extends into a corona above (at least the inner por¬ 
tion of) the thin accretion disk. We assume that there is 
significant flare activity produced through magnetic re¬ 
connection within this corona (e.g. Kuijpers 1995). Such 
releases of magnetic energy are presumably concentrated 
near the inner boundary of the disk where the densi¬ 
ties and field strengths withi n the disk are greatest (e.g. 
Ghakrabarti & D’Silva 1994 ), and thus we will consider 
parameters relevant to this region. 

Magnetic reconnection is now believed to be an es¬ 
sential component of the solar flare process (see Lakhina 
200C for a recent review). Bi-directional electron beams 
(i.e., propagating both towards and away from the Sun) 
are now understood to be a signature of the magnetic re¬ 
connection process as revealed throu gh solar radio obser¬ 
vations (e.g. Robinson & Benz 2000). Here we investigate 
conditions under which the “chromospheric evaporation” 
induced by downgoing electron beams (e.g. Fisher 1987; 

could operate in ADAF type accre- 


Bastian et al. 

tion disks giving rise to rapid variability or microflares. 
In analogy with the solar case, we anticipate that a large 
number of small flares will be present simultaneously (e.g. 
Parker 1988). The evaporation of thin disk material can 
erode the inner portion of the disk, thus changing rt, and 
thereby possibly inducing additional X-ray variability. 


3. Evaporation from accretion disks 

3.1. Conditions appropriate to low states for XRBs 


To illustrate this situation, we adopt simple standard a 
disk models for the outer port ion of the accretion disk 
(e.g. Shakura & Sunyaev 1973 , hereafter SS), but trun¬ 
cate them at rt, where the two-temperature advection 
dominated portion of the flow is assumed to take over. 
In order to incorporate the explosive gasb ag (E G) model 
for chromospheric evaporation (e.g. Fisher 1987 ) we need 
values for the density rich and the scale height, zo in that 
region. At this distance from the BH, the thin disk will 
be characterized by gas pressure dominance over radia¬ 
tion pressure and electron scattering opacity larger than 
free-free or other opacities, so (SS), 


h = 2.3x 


no = 9.2 X 


- 7 / 10 ^- 7 / 10 ^ 2 / 5 ^- 


33/20cm-3. 


( 1 ) 

( 2 ) 


Here, h is the disk half-thickness, a the viscosity param¬ 
eter, m = Mbk/Mq, fn = M/{2.3 x lO“®M 0 yr“^m), 



BH 


ADAF region 



Fig. 1. A schematic of an accretion flow in the low-state, 
where the transition radius, rt, is far from the black hole 
(BH). The inner, two-temperature plasma extends out to 
rt, with a thin disk and a corona present beyond that 
point. Flares within the corona produce electron beams, 
some of which impinge upon the cool disk and evaporate 
material from it into the corona. 


r = R/ (2GMbhc“^), no is the number density of the disk 
in the equatorial plane, and R is the physical distance from 
the BH, s o the coefficients differ from those in Shakura & 
Sunyaev (1973), which uses different definitions of rh and 
r. We have dropped the relativistic correction terms which 
are negligible at the distances considered here. 


These standard disk models do not provide precise val¬ 
ues for the density, nch and the scale height, zq, in the pho- 
tospheric/chromospheric regions of the disk. Nonetheless, 
an estimate for the zq in the photosphere, where the den¬ 
sity falls exponentially as n ex exp(— z/zq), is given by 
Zo = r^kT{zi)/(GMnripZi), where rup is the proton mass 
and zi is given by the location where the total optical 
depth equals unity (SS). In that Zi is only slightly greater 
than Zo and the chromosphere begins at a height of at 
least 2 zi, we adopt values of rich = 1-0 x 10 “^ y no and 
Zo = 0.3 rj h, where we expect y, and rj to be within an 
order of magnitude of unity. 
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The microflares are likely to be concentrated in the 
transition region, and we use parameters relevant for low- 
states of ADAF models. This implies that m < 1.3a^ 
(e.g. Esin et al. 1997). Most current estimates for disk 


where rup is the proton mass, and Zq is the initial thickness 


o A in5 -1/10 9/10 . 1/5 21/20 

zq = 3.4 X 11 ) 77 a ^ m /i rmn cm. 


.01' 100 


Uch = 7.3 X 10^®x a 1 


-7/10 -7/10 . 2/5 -33/20 


''10 


m 


01' 100 


cm 


( 3 ) 

( 4 ) 


gasbag model (e.g. Fisher 1987). Using the standard equa- 


( mpZ'^ Y/3 

\ ttQ ) 


ttQ 


(6) 


viscosity involve some form of the magnetorotational in¬ 
stability and tend to produce a < 0.1 (e.g. Brandenburg 
199S; Hawley et al. |1996), and we thus assume typical 


values of a,i = a/0.1 ~ 1, rh.oi = m/0.01 ~ 1. For these 
low states, rt ~ 10^ — 10^, so we parameterize by taking 
T’loo = rt/100. We shall first focus on LMXBs, for which 
a typical BH mass of IOMq can be assumed, and thus 
use mio = Mbh/(10A'/q) ~ 1. Then we obtain the fol¬ 
lowing estimates for the scale height and density in the 
chromospheric regions 


of the evaporating plasma. As shown by Fisher (1987), if 
this EG model is to work at all, Zq must differ from zq 
only by a logarithmic factor not much larger than unity. 
For the circumstances considered here, conduction is not 
important until times exceeding ~ lOr, and so the time r 
represents the characteristic time over which fluctuations 
are produced by this process. Of course, other processes 
which we are not considering here may well play signihcant 
roles and could produce flares by different mechanisms on 
different timescales. Identifying Zq = (^zq with (/ > 1, we 
have the following approximate expression 


= 2.3 X s. 


( 7 ) 


During a reconnection event, electron beams form in 
the current sheets; some of these will be directed down¬ 
wards into the disk atmosphere and heat it. Radiative 
losses must be considered in this situation. In this region of 
the accretion disk the standard photospheric temperature 
is given by 

r, = 1.4 X K. (5) 

The radiative loss function, A(r), is maximized at tem¬ 
peratures near 2 x lO^K (Raymond et al. 1976| ), which is 
marginally above the value just found for Tg. Thus it is 
likely that such temperatures will be found in the tran¬ 
sition region between the photosphere and the corona of 
the accretion disk for rioo > 1, and this maximum loss 
rate can be achieved. In order for this evaporation mecha¬ 
nism to function, the minimum heating rate per particle, 
Q, must at least balance the maximum cooling rate, and 
is foun d to b e Q = A(r)uch — 7 x 10“^^erg s“^cm“^ rich 
(Fisher 1987), or 


Therefore, for XRBs, timescales of less than a millisecond 
are possible for the explosive evaporation from the cool 
disk, even in the low state, where rt is at least 100. 

The minimum electron flux needed for balancing the 
radiative losses can be approximated by 


F cx 




Q = 5.1 X 10 aoi''®riOQ^''^°erg s ^ 

If the photospheric temperature exceeds 2 x lO^K then the 
cooling rate drops rapidly to about Q = 1.5 x 10“^^nch 
by Tg = 5 X lO^K, and then much more slowly until Tg ~ 
3x 10®K (Raymond et al. 1976). Under these cicumstances 
heating rates several times lower than that given above 
can suffice. 

The energy deposited by the electron beam(s) heats 
the chromosphere, part of which evaporates high into the 
corona. If the other transport processes (e.g. conduction) 
are not immediately significant, then the evaporation can 
occur in an explosive manner and can be described by the 


tions for mass, momentum and energy conservation, one 
can derive the time prohles of the plasma parameters in 
this gasbag. Of key importance is the timescale for build 
up of maximum pressure and consequent rapid evap ora- 
tion of the plasma, which is found to be (Fisher 1987) 


where Nc = 9.15 x 10^®cm ^(Ecpo)^ (McClymont & 
Canfield 1986| ; Fisher et al. 1985). Here, S is the power- 
law index for the non-thermal distribution of beam elec¬ 
trons (in the sense that f{E)dE oc E~^dE), and the solar 
flare obs ervati ons indicate that 5 ~ 4 — 5 (McClymont & 
Canfield 1986| ). The total column density of beam elec¬ 
trons, Nc, is determined from collisional processes and 
depends upon the cut-off energy to which they can pen¬ 
etrate, which is typically 5-20 keV, so we use Ec^w = 
Ac/lOkeV ^ 1. This electron energy flux E implies densi¬ 
ties in the beam several orders of magnitude lower than the 
ambient density, rich, so the standard acceleration mecha¬ 
nisms invoked in solar flares are adequate (e.g. Bastian et 
al. 1998 ). Explosive evaporation occurs for fluxes greater 
than this minimum value, and this EG process is not rel¬ 
evant if the fluxes are less than this. 

To obtain the total power emitted in a microflare 
we must estimate the area over which this flux emerges. 
Again, in analogy with the sun, we approximate the size, 
R, of the flaring region to be of the order of the scale 
height, Zq. However, it is certainly possible that the flaring 
loop could extend into the corona far enough so that the 
appropriate size is actually that of the coronal scale height, 
which would be several times larger (e.g. McClymont & 
Canfield 19861 ). Therefore, since nearly all of this power 
will be quickly radiated in EUV and X-ray photons, we 
can obtain an estimate for the minimum total power in 
one of these microflares as, L/ = ttR^F; taking 5 = 4 we 
have, 


Lf ~ b.lxlQ^^rfxoi.i 


-9/10 11/10 . -4/5 9/20 

'"-10 '"o.oi 'lOO 


-(f) 


erg s 


The power in a single microflare Lf\s more than six orders 
of magnitude smaller than the total power in the low state, 
using our nominal values of 77 cx x — (R/zq) ~ 1, so the 
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sobriquet microflare is quite appropriate. For these low 
state conditions, (r^ > 100), each of these nano/micro¬ 
flares covers an area very small in comparison with the 
active portion of the disk (~ 10“®). It is essentially im¬ 
possible to properly estimate the fractional area of the 
relevant part of the accretion disk that will be subject 
to this process at any particular time. However, it is cer¬ 
tainly plausible to expect that no more than 10% of the 
disk should be so involved, and thus it is unlikely that 
> 10® of these explosive events can be going on simulta¬ 
neously. Even with this optimistic assumption, the inte¬ 
grated output of all these microflares is likely to remain 
small in comparison with the basal luminosity. Only if 
most of the variables in Eq. (8) exceed their nominal val¬ 
ues could these microflares make a significant contribution 
to the total power. Probably the least certain factor in Eq. 
(8) is the ratio R/zq, and it is not implausible that this 
geometrical factor could approach 10; in that case, the in¬ 
dividual microflare luminosities rise to above 10“® of the 
observed total low state power. Only in this optimistic sit¬ 
uation is it possible that the integrated flux from all such 
flares could amount to a substantial fraction of the total 
emission. However, for large values of i?, we expect there 
to be fewer of these “miniflares”, in that each now covers 
a larger fraction of the disk area. Hence we do not expect 
the integrated micro/miniflare power to rise as steeply as 

3.2. High/soft state conditions for XRBs 

Here we consider the situation where the accretion rate is 
high enough so that the standard accretion flow penetrates 
much closer to the BH. Then the relevant equations are 
those appropriate for the inner portion of accretion disks, 
which are dominated by radiation pressure (SS). We note 
that for a high state, m ^ 0.1 and r* ~ 10 are typical 
parameters; furthermore, the relativistic corrections are 
no longer negligible, since we are dealing with regions close 
to the BH. With the reparameterizations m,i = m/0.1 
and rio = r/10, defining B = [1 — (0.3/rio)^^^], and still 
adopting zq = O.Sry/i and rich = 1.0 x 10“^xno, we find: 

zq = 7.4 X 10'^ r)miQrn,iB cm, (9) 

rich = 4.4 X 10^®XQ:”)^m/g^m“)^rJg^H“^ cm“®. (10) 

Under these circumstances the photospheric tempera¬ 
ture is given by 

Ts = 6.6 X 10®aY^m”^/®mY'^r“^''^H4/® K, (11) 

or Ts ~ 3.5 X 10®K for rig = 1, and all other quantities 
taking on their new nominal values. This somewhat ex¬ 
ceeds the temperature at which maximum cooling occurs. 
It remains reasonable to utilize the value of A correspond¬ 
ing to this higher temperature to compute the heating 
rate needed to compensate for cooling in the optically thin 
regime; for this temperature, A is roughly half that of the 
peak. 


Proceeding as in the previous subsection we then ob¬ 
tain, 

Q ~ 1.5 X 10“^ xa^^m/g^m^^rjg^i?“^erg s“^, (12) 

T = s, (13) 

and 

Lf = 7 X 10^®?7^xa“i^mio erg s“^ (14) 

In that the total power for an XRB in the high state is ~ 
10®^erg s“^, we see that each of these micro (nano?) flares 
contributes very little to the total power output. Of course, 
as for the low state, many such small scale flares may 
contribute simultaneously; however, in that the relevant 
disk area under the intense coronal region is proportional 
to and is thus smaller than in the low state case, the 
integrated contribution is very likely to be negligible. 


3.3. Results for AGN 


Active galactic nuclei may also be described within the 
ambit of ADAF models, particularly those with small out- 
put , and therefore proba bly lo w efficiency (e.g. Lasota et 
al. 1996 ; Narayan et al. 1996| ) so essentially everything 
developed above can still be applied. 

The photospheric temperatures will be significantly 
lower, since they decline with increasing black hole mass, 
but coronae and transition regions should still be present, 
and microflares can still evaporate disk material into the 
corona. Appropriate values for all dimensionless param¬ 
eters are nearly the same as in the two preceding sub¬ 
sections. Therefore the timescales vary as M®/® for the 
low state and as M for the high state. Roughly, then, 
these microflares operate on 3-10,000 s timescales (for 
Mbh ~ 10^ - IO^Mq). 

Since A/ oc for the low state and Lf or. M for 

the high state, each of these microflares could contribute 
between ~ 10®"^ — 10®® erg s“®, while the total powers are 
~ 10'^® erg s“® for low-state emission around a IO^Mq 
BH up through ^ 10"'^® erg s“® for high-state emission 
near a 10® solar mass black hole. Thus, even considering 
the possiblity that many of these microflares can erupt 
simultaneously, their contribution to the total power still 
should be very small. 


4. Conclusions 

The mechanism of explosive chromospheric evaporation 
seems to explain some aspects of solar X-ray variability, 
as well as blue-shifted soft X-ray lines from the Sun. Since 
conditions in accretion disk atmospheres probably have 
many similarities to the solar case, we have explored this 
mechanism in the context of accretion disk models for 
XRBs and AGN. In particular, we have focused on the 
transition region in ADAF models as the most likely place 
for these flares and microflares to be generated. We have 
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estimated microflare timescales and powers based on the 
exploding gas bag model for both low and high states for 
both XRBs and AGN. 

We note that we have not proven that this mecha¬ 
nism is necessarily relevant; we have merely shown that 
it is plausible that it could operate in the conditions of 
accretion disk atmospheres, as it probably does in the 
Sun. Furthermore, it is very likely that other flare mech¬ 
anisms dominate the observed variability. For XRBs, the 
EG timescales range from sub-ms to seconds, while for 
AGN, timescales from seconds to hours are expected, with 
most of the relevant emission emerging up to ~ 10 keV. 
Whereas these characteristic times are in the range of 
the fastest observed variability for both classes of sources, 
the powers in individual microflares are extremely small. 
However, we do expect a large number of these flares could 
exist essentially simultaneously, particularly in the low 
states. Even so, the integrated emission from these mi¬ 
croflares is almost certainly undetectable with the current 
generation of X-ray telescopes. 
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